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A defect approach coupled with matched asymptotic ex- 
pansions is used to derive a new set of boundary layer equa- 
tions. This method ensures a smooth matching of the bound- 
ary layer with the inviscid solution. These equations are 
solved to calculate boundary layers over hypersonic blunt bod- 
ies, involving the entropy gradient effect. Systematic compar- 
isons are made for both axi symmetric and plane flows in sev- 
eral cases with different Mach and Reynolds numbers. After 
a brief survey of the entropy layer characteristics, the defect 
boundary layer results are compared with standard boundary 
layer and full Navier- Stokes solutions. The entropy gradient 
effects are found to be more important in the axisymmetric 
case than in the plane one. The wall temperature has a great 
influence on the results through the displacement effect. Good 
predictions can be obtained with the defect approach over a 
cold wall in the nose region, with a flrst order solution. How- 
ever, the defect approach gives less accurate results far from 
the nose on axisymmetric bodies because of the thinning of 
the entropy layer. 


Introduction 

A blunt body in hypersonic flow is preceded by a bow 
shock wave, detached in front of the nose. This strong curved 
shock wave induces an entropy gradient in the shock layer. 
For an inviscid flow, the entropy gradient is related to the 
vorticity through Crocco equation ; 

curl V A V = — grad Ht grad S 

Therefore, velocity and temperature gradients also exist in 
an inviscid shock layer. The standard boundary layer theory 
of Prandtl cannot take into account normal gradients out- 
side of the boundary layer. Van Dyke proposed an enlarged 
theory called higher-order boundary layer theory based on 
matched asymptotic expansions for high Reynolds numbers 
[11, 12]. Two expansions corresponding to different approxi- 
mations of the Navier-Stokes solutions are built. One of them 
called outer expemsion is valid far from the wall, where the 
viscous effects are negligible. The other one called inner ex- 
pansion describes the boundary layer where the viscous effects 
are dominating. The Prandtl boundary layer equations then 
represent the first term of an expansion in powers of a small 
parameter. The external flow normal gradients are accounted 
for in the second order term, which is a small perturbation 
of the first order solution. Several other second order effects, 
like the wall curvature, the displacement or the rarefied gas 
effects are brought into evidence. The main advantage of this 
systematic method is to give not only the equations but zdso 
the matching conditions between the different zones. 


For hypersonic reentry flows, the Reynolds number is of- 
ten moderate at high altitudes, because of the low density of 
air. The boundary layers are thus thick and can be of the 
same order of magnitude as the entropy layer, and the invis- 
cid flow quantities can undergo important variations between 
the wall and the edge of the boundary layer. Because of the 
hypothesis of Reynolds number tending towards infinity, the 
boundary layer is assumed to be very thin in Van Dyke’s the- 
ory and the inviscid flow gradients are represented only by 
their wall value. So the second order expansion cannot ensure 
a good matching of the boundary layer with the inviscid flow 
if the inviscid profiles are not linear, and the influence of the 
external vorticity on the skin friction and the wall heat flux 
is not correctly estimated. 

Defect approach 


Decomposition 

To ensure a smooth matching at any order whatever the ex- 
ternal flow, a defect approach has been used, coupled with 
asymptotic expansions [3j. In the boundary layer region, the 
variables are no longer the physical variables, but the dif- 
ference of them with the external solution (Le Balleur [9]). 
We consider a steady two-dimensional flow of ideal gas. The 
variables /?, u, u, p and T stand for the density, tangential 
and normal component of the velocity, pressure and temper- 
ature. The equations are written in a system of orthogonal 
curvilinear coordinates where ( represents the curvilin- 

ear abscissa along the body and rj is the distance to the wall. 
All the variables are made dimensionless by referencing them 
to the upstream values and Uqo, the nose radius Rq and 
To = Ul^lCp. So we write : 

P = Pe Pd 


U = Ue -h 

V = -f U£> - 0) 

P = Ps + Pi? 

T ^Tb + Td 

where the subscript E stands for the external variables and 
the defect variables are labelled D. The term ve(^, 0) ha^ 
been added to keep the condition up(^,0) = 0 at the wall 
whatever the value of vg. 

Expansions are then written using the same small param- 
eter e as Van Dyke : 


c 


1 


Re = 


PaaUaaRa 

KTo) 


The external functions depends on the coordinates The 

outer expansions read : 

+ SU2U,V) + ••• 

= Vi{i,Ti)/+ cVliCv) + ••• 

Pe{(,v) = + sP2{Cv) + 

peU,v) = RiU,^) + + ■■■ 

Te{(,v) = Pi{(>v) + + ■■■ 

In the inner region, a stretched normal coordinate ^ = tj/c is 
used for the defect variables : 

^dU, ri) =ui{(,rj) +euj{(,^) +■■• 

vdU.»?) = ev,(4,^) + e’v3(^,77) + 

pdU,!) = +^Pi{(,n) + ■•• 

PD{i,v) = Piii,v) +spj{Ln) + ■■• 

PdU,v) = + ■■■ 

The expansion for v must be shifted to avoid the degenercury 
of the continuity equation. These expansions are then brought 
into the Navier-Stokes equations, and terms of like power of 
€ are equated. 


- state : 

Pi = — [Pi^i (^1 Pi)^i ’ 

The symbol r represents the distance from the wall to the 
symmetry axis, with ; = 0 for plane and j = 1 for axi sym- 
metric bodies. As in Prandtl equations, the wall curvature 
appear in the first-order equations only through the trans- 
verse curvature radius in the continuity equation. The second- 
order equations are small- perturbations of the above ones plus 
source terms due to curvature effects, like in Van Dyke theory. 

Matching conditions 

Each expansion must satisfy the boundary conditions corre- 
sponding to its own validity domain. The upstream conditions 
are to be applied to the outer expansion and the wall condi- 
tions to the inner one. The missing conditions are obtained 
by matching the inner and outer expansions. At the edge of 
the boundary layer, we can write : 

u — ► us 

U -♦ V£; 

P ^ PE 
P ^ PE 

T ^Te 

and so for the defect variables : 


variables are null and the equa- 
tions for the outer flow arc exactly the same as for Van Dyke’s 
theory, i.e. Euler equations. Concerning the inner region, one 
must first bring the above expansions into the Navier- Stokes 
equations, then subs tract the external equations, and at last 
equate same powers of e. For practical convenience, the inner 
equations can then be rewritten in outer coordinates, using rf 
instead of 77 , and replacing Vi and by : 


Equations 

In the outer region, the < 


Then the following first-order equations are obtained : 
- continuity : 




drj i 


dvi 


dr\ 




-momentum : 

dU 

[Ri -h P\){Ui -H ^1)“^ + \pi^i + {^1 + 


+ Pi)(K + Vi) 


dui 


5pi ^ 

Re St; dr} ] 


- r^-momentum : 


0 = 


dp\ 

ft? 


- energy : 

3t ST 

{Rr+Pl){U2+U,)-^ + \P2U2 + {Rx + 

+(Ar + /h)(Vx + + (C;, + 

^St; VPr Re Sr;/ Re V St?/ 


Uij 0 

vd ^ ve(^,0) 
Pd — ^ 0 
Pd 0 

Td->0 


Thus at first order : 


lim ui = 0 

^J-*oo * 


V^iU,0) = 0 


lim Pi = 0 
n— 00 ^ 

lim ti = 0 


lim Pi = 0 


The conditions on p, p and T are not in dependant since they 
are linked through the state equation. The condition on v is 
not a boundary condition for the inner expansion but it gives 
the wall condition for the outer flow. 

The wall conditions for the inner flow are : 


u = I/i Ui -f- €{ 1/2 + Tij) = 0 
V = CUi H- £^V2 = 0 

T = Tj -h ti + c{T2 + tj) = T^ 


wi(^O) 

wi(^,0) 

tiU,0) 


0 


hence : 



Discussion 

Thanks to the small perturbation approach, the calculations 
of external flow and boundary layer are uncoupled and can 
be performed separately provided that a specified sequence 
is respected. First order external problem must be solved 
first, then first order internal, second order external, and so 
on. The defect boundary layer equations are parabolic and 
C8U1 be solved by space marching at a very low cost, like the 
standard Prandtl equations. 

The conditions at the edge of the boundary layer ensure 
a smooth merging of the boundary layer into the inviscid flow 
whatever the inviscid profiles. From a theoretical point of 
view, it can be shown that the defect expansions are consistent 
with Van Dyke’s ones by the fact that at a given order they 
differ only by terms which are higher-order in Van Dyke’s 
theory. 

Using the above conditions, the first order 77-momentum 
equation reduces to 

Pi = 0 

So, the pressure in the first-order boundary layer is every- 
where equal to the local inviscid flow pressure, instead of its 
wall veilue like in Van Dyke’s theory. 

Applications 

To experiment the defect approach, several cases have been 
selected for a blunt body in a hypersonic flow of ideal gas. 
The general shape of the body is a plane or axisymmetric 
hyperboloid, defined by the nose radius and the angle of the 
asymptotes, at zero degree incidence. The numerical data are 
given by Shinn, Moss and Simmonds [10] for a hyperboloid 
equivalent to the windward symmetry line of the U.S. space 
shuttle. Two points of the reentry trajectory of the STS-2 
flight are considered here : 


Reentry trajectory - Flight STS-2 

Mach Mqo 

26.6 

23.4 

time (s) 

250 

650 

altitude (km) 

85.74 

71.29 

nose radius Rq {m) 

1.322 

1.253 

asymptotes half- angle (°) 

41.7 

40.2 

pressure poo (P*) 

0.3634 

4.0165 

temperature To© (K) 

199 

205 

velocity Uao (ni/s) 

7530 

6730 

density p„ (kg/m^) 

6,35 10'® 

6.80 10-® 

reference temperature To (K) 

56321 

44900 

Reynolds number Re = 

183.55 

1865.65 

small parameter e = Rc”*^* 

0.074 

0.023 

Reynolds number Rcoo = 

4792 

42374 

Moo 




The Prandtl number is assumed to be constant and equal 
to 0.725. The ratio of specific heats 7 is 1.4. The wall tem- 
perature is fixed and equal to 1500 K. The viscosity law is 

Sutherland’s. No comparison with experimental data is pos- 
sible since the real gas effects are not yet included. So Navier- 
Stokes solutions [8] have been taken as reference, to compare 


the two Euler + boundary layer methods. Euler calculations 
are made with a code from ONERA [14]. Standard bound- 
ary layer solutions are obtained using a program developed 
in DERAT [2j. Only first-order boundary layer are presented 
here since second-order outer flow solutions are not yet avail- 
able, Several second-order calculations using Van Dyke’s the- 
ory have been made on a hypersonic blunt body [1, 4, 5, 6, 7]. 

Axisymmetric hyperboloid 

Past a hyperboloid, the shock wave curvature decrease fastly 
and the entropy field tends to be uniform, except for the 
streamlines near the wall, which crossed the strongly curved 
shock wave at the nose. In this case, the entropy layer is 
characterized by a non- zero normal gradient at the wall and 
a decreasing thickness towards the rear, since the mass- flow 
is constant in the entropy layer and the circumference of the 
body increases (fig. 1 left). The entropy values at the wall 
and at the edge of the entropy layer remain constant because 
the wall is a streamline and outside of the entropy layer the 
flow is isentropic. So the normal entropy gradient at the wall 
deeply increases downstream. The shock layer is thinner than 
in the plane case. Far from the nose, the flow is similar to 
a flow past a sharp cone except in the entropy layer, whose 
aspect is quite similar to a viscous boundary layer (fig. 2). 

Boundary layer profiles are displayed on figures 4 to 7 for 
the Mach 23.4 case. Longitudinal velocity profiles are plotted 
on figure 4 at a distance of nine nose radius from the stagna- 
tion point. One can see on this figure the important velocity 
gradient at the wall in the inviscid flow. This gradient dimin- 
ishes distinctly between the wall and the boundary layer edge. 
So even with a second-order expansion, Van Dyke’s method 
could not give a good matching, since it considers only the 
wall value of the gradient. In this case, it would widely over- 
estimate the skin friction (Adams [l]). Due to the very lov: 
wall temperature compared to the inviscid flow one, the dis- 
placement effect is quasi-null and the Navier- Stokes solution 
recasts exactly the inviscid profile in the outer region. In this 
case, the agreement is quite good with the first-order defect 
boundary layer. A composite profile has been plotted also, 
using the additive composite expansion (Van Dyke [13]) con- 
structed with the first order inner and outer expansions. It 
gives good results for the longitudinal velocity, slightly differ- 
ent of the defect ones. 

The corresponding profiles for the temperature are shown 
on figure 5. The defect profile is in rather good agreement 
with the Navier-Stokes solution, but in this case the compos- 
ite expansion written with Van Dyke’s first order solutions 
gives very bad results and does not improve the inner solu- 
tion. This is due to the negative slope at the wall for the 
inviscid temperature. Figures 6 and 7 show the velocity and 
temperature profiles at twenty- one nose radius. The growing 
boundary layer has overlapped a larger part of the entropy 
layer. Because of the constant total enthalpy, the positive 
velocity gradient at the wall induces a negative temperature 
gradient. In spite of this, the wall heat fliix is increased by 
the vorticity, as well as the skin friction, as can be seen on the 
figures 8 and 9. But the increase is far more important for the 
wall friction than for the flux. The defect approach underes- 
timates slightly these quantities but gives better predictions 
than the standard boundary layer. 





02468 02468 


axisymmetric plane 

Figure 1: Entropy level in the shock layer - Mach 23.4 



Figure 3: Entropy profile - plane case Mach 23.4, { — 9 


Figure 2: Entropy profile - axisymmetric case 
Mach 23-4, ^ = 9 










skin friction 






Figure 10; Longitudinal velocity profiles 
Mach 23.4, = 15000 K, f = 4 


Figures 10 and 11 show the velocity and temperature pro- 
files with an arbitrary temperature of ten times the temper- 
ature of the preceding case. The displacement effect is then 
far more important and it is obvious on these figures that the 
Navier- Stokes solution is shifted from the Euler solution in 
the outer zone. So the first-order boundary layer methods 
give poor results and a second-order calculation seems to be 
necessary. 

The velocity and temperature profiles at nine nose radius 
abscissa in the Mach 26.6 case are presented on figures 12 
and 13. Because of the lower density, the Reynolds number 
is small eind the boundary layer is about twice as thick as in 
the Mach 23.4 case. So a large part of the entropy layer is 
overlapped by the boundary layer. The inviscid velocity and 
temperature gradients at the edge of the boundary layer are 
far weaker than their wall values. Due to the high value of the 
expansion parameter e, the second order effects are more im- 
portant and a slight displacement effect is visible between the 
Euler and Navier- Stokes profiles outside the boundary layer. 
The agreement between the Navier- Stokes and defect profiles 
is rather good, but the shear at the wall is a bit too high 
for the later one. Note that because of the negative inviscid 
temperature gradient at the wall, the Van Dyke’s composite 
expansion gives again poor results on the temperature profile. 

Figures 14 and 15 show the same quantities at twenty- 
one nose radius from the nose. The entropy layer is now com- 
pletely included into the boundary layer, and the gradients 
in the entropy layer become higher than those of the viscous 
boundary layer. So the hypothesis of neglecting the viscous 
effects in the external flow does not hold any longer and the 
defect boundary layer probably gives overestimated values for 



Figure 11: Temperature profiles 
Mach 23.4, = 15000 K, ^ = 4 


the slope at the wadi of the velocity profile. But no Navier- 
Stokes solution is yet available on such a large domain. 

The corresponding skin friction and wall heat flux are 
shown on figures 16 and 17. As forecast from the velocity 
profiles, the defect approach improves greatly the standard 
boundary layer result, but widely overestimates the skin fric- 
tion on the rear of the body. The predictions concerning the 
wall heat flux seem to be more reliable. 

Plane hyperbola 

Let us now consider a plane hyperbola in the same conditions 
of hypersonic flows. On figure 1- right are displayed the en- 
tropy levels in the inviscid shock layer. The main difference 
with the axisymmetric case is that now the entropy gradient 
is null at the wall (Veui Dyke [12]). Figure 3 shows entropy 
profile across the shock layer. The entropy gradient layer is 
thus located at a short distance above the wall. So the ve- 
locity and temperature gradients in the inviscid flow are null 
at the wall as well, and their influence will be significant only 
with a very thick boundary layer. Moreover, far downstream, 
the flow can be assimilated to a parallel flow and the entropy 
layer’s thickness remains constant whereas in the axisymmet- 
ric case the entropy layer gets thinner towards the rear of 
the body. Thus the entropy gradient remains bounded. Since 
it is null at the wall, its influence on the skin friction and the 
heat flux will now be far less important. 

On figures 18 and 19 are plotted the velocity and temper- 
ature profiles on the Mach 23.4 hyperbola at nine nose radius 
abscissa. The inviscid gradients are hardly visible outside the 
boundary layer and all the methods give the same results. 

When the Reynolds number is lower, the matching of the 




Figure 14: Longitudinal velocity profiles Figure 15: Temperature profiles 

Mach 26.6, T, = 1500 K, ^ = 21 Mach 26.6, T, = 1500 K, f = 21 
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skin friction 



Figure 16: Skin friction on the hyperboloid - Mach 26.6, Ti, = 1500 K 



Figure 17: Wall heat flux on the hyperboloid Mach 26.6, = 1500 K 
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Figure 18: Longitudinal velocity profiles 
Mach 23.4, = 1500 K, ^ = 9 



Figure 20: Longitudinal velocity profiles 
Mach 26.6, = 1500 K, ^ = 9 



Figure 19: Temperature profiles 
Mach 23.4, T* = 1500 K, { = 9 



Figure 21: Temperature profiles 
Mach 26.6, = 1500 K, { = 9 





boundary layer with the in viscid flow takes place in the gra- 
dient region, as can be seen on the flgures 20 and 21 for the 
case Mach 26.6 . The defect method gives a good matching 
and a correct agreement with Navier-Stokes solutions, but the 
two boundary layer methods give similar results near the wall. 
Thus no significant difference is visible on the skin friction and 
the wall heat flux shown on figures 22 and 23. 

Conclusion 

Several boundary layer calculations have been performed on 
various hypersonic bodies, including plane or axisymmetric 
shapes. The behaviour of the solutions far from the stag- 
nation point has been particularly investigated. The different 

cases presented here showed both the interest and the limits of 
boundary layer methods to compute hypersonic flows. More- 
over, the importance of taJeing into account the second or- 
der effects when calculating boundary layers at low Reynolds 
numbers has been brought into evidence. The most important 
of them are the entropy gradient effect and the displacement 
effect in the considered cases. They can deeply modify the 
wall quantities_su^ the skin friction or^t he wall heat flnx, 
which are essential to predict the total drag of the vehicle and 
to design the thermal protection. 

Using the matched asymptotic expansions technique, the 
defect approach allows us to improve the results of the stan- 
dard higher- order boundary layer theory of Van Dyke, for a 
similar cost. Particularly, it ensures a smooth matching of 
the viscous and in viscid flows, even when the in viscid profiles 
vary significantly through the boundary layer. When the wail 
temperature is low and thus the displacement effect is negli- 
gible, first-order defect calculations can give good results and 
reproduce Navier-Stokes solutions with a reasonable accuracy 
at a lower cost, as long as the entropy layer is not too thin 
compared to the viscous boundary layer. But it gives less ac- 
curate results on axisymmetric hyperboloids at low Reynolds 
number far from the stagnation point, wh_en the in viscid flow 
normal gradients are higher than those of the boundary layer. 
The inviscid flow concept seems to be invalid then. 
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